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^ The recent discovery of J0737-3039A 1 & B 2 -two pulsars in a highly relativis- 

o : 

O ; 

the elusive physics of pulsar radio emission. The system contains a rapidly 
rotating pulsar with a spin period of 22.7 ms and a slow companion with 
a spin period of 2.77 s, hereafter referred to as 'A and 'B, respectively. A 



tic orbit around one another - offers an unprecedented opportunity to study 



unique property of the system is that the pulsed radio flux from B increases 
systematically by almost two orders-of-magnitude during two short portions 
of each orbit 2 . Here, we describe a geometrical model of the system that 
simultaneously explains the intensity variations of B and makes definitive 
and testable predictions for the future evolution of the emission properties 
of both stars. Our model assumes that B's pulsed radio flux increases when 
illuminated by emission from A. This model provides constraints on the 
spin axis orientation and emission geometry of A and predicts that its pulse 
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profile will evolve considerably over the next several years due to geodetic 
precession until it disappears entirely in 15-20 years. 

The double-pulsar system is extremely compact (orbital period, P orb ps 2.45h), 
mildly eccentric (e « 0.088), and highly inclined ( % > 87°). The orbital light curve 
shows two regions that last for about 12 min (30° of orbital phase) and are separated 
by approximately 74° where the pulsed radio flux of B increases dramatically 2 . The 
average pulse profile of A is composed of two peaks of width 20° — 30° separated by 
about 150° in pulse phase 1 . 

We report on a model that explains the pulsed intensity variations seen in the orbital 
light curve of B using two simple assumptions. First, we assume that the emission 
geometry of A is given by the traditional hollow cone model of Lyne & Manchester 3 
(see figure 1), and second, that B appears bright only when the emission beam of 
A illuminates it (i.e. when the cone of emission from A intersects B's position in the 
orbital plane). The key feature of the second assumption, which we call the "stimulated 
emission hypothesis", is that the particles and/or electromagnetic radiation from A 
which stimulate the radio emission from B are directed along the same emission cone 
as A's observed radio radiation. These two assumptions enable one to derive a set 
of non-linear equations (see the caption to figure 1) that determines how the conal 
emission from A intersects with our line-of-sight, how it projects onto the orbital plane, 
and when it illuminates B. These five equations contain a total of 11 variables, six of 
which are obtained from observations, and which therefore allow us to solve for the five 
unknowns. The six measured parameters are the orbital inclination angle, i, the angle 
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between the outer edges of the double peaked pulse profile, P Q , the angle between the 
inner edges of the peaks in the pulse profile, Pi, the angle between the outer edges of the 
double peaked orbital light curve, Q , the angle between the inner edges of the orbital 
light curve, Oi, and the midpoint of the orbital light curve, O m . The angles obtained 
by solving the five constraint equations are A and <p which are the polar angles that 
determine the orientation of the spin axis with respect to the orbital angular momentum 
(see figure 1), and p, a, and 5, which determine the emission geometry (see figure 2). 

Table 1 lists the adopted values of the measured parameters i, Pi, P , O h O a , and 
O m . These values were obtained from 427, 820, 1400, and 2200 MHz data taken with 
the Green Bank Telescope 4 . These observations, which cover a factor of 5 in radio 
frequency, indicate that the pulse profile of A and the orbital light curve of B are both 
effectively independent of observing frequency for our purposes. Given the measured 
parameters, we solved the constraint equations and found two solution sets for the five 
unknown angles, which are also listed in Table 1. These two solutions exist due to the 
uncertainties in the measured parameters. A video showing the relative motions of the 
pulsars, the intersection of A's conal emission with the orbital plane for solution 1, and 
the regions of enhanced emission from B is available on-line 5 . 

In our model, the orientation of A's spin axis plays a major role in determining 
the shape of both the pulse profile of A and the orbital light curve of B. Geodetic 
precession 6 ' 7 , caused by the curvature of space-time around the pulsars, will cause the 
spin axis to precess about the orbital angular momentum vector, J, at a rate of 4.77° per 
year. Hence, the angle <p will increase at this rate while the four other inferred angles 
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will remain constant. Consequently, we can predict the future evolution of both the 
orbital light curve and the pulse profile of A. Given the inferred values of A, 0, p, a and 
S together with a value of at a specified time, the five constraint equations determine 
the pulse profile parameters, Pi, and P Q , and the orbital light curve parameters O a , 
Oi, O m . Both Oi and Q will remain constant while O m , which is equal to 0, will 
increase at a rate of 4.77° per year. This effect will be easily measurable over the 
next 1-2 years. The geodetic precession will also cause highly significant pulse profile 
changes over timescales of a few years (see figure 3). Note that similar, although smaller 
amplitude, precession-induced profile changes have already been observed in the binary 
system PSR B1913+16 8 ' 9 . For solution 1, P and P { will increase by 42° ± 16° and 
31° ± 8° in 1 year , respectively. For solution 2, P Q and Pj will increase by 96° ± 13° and 
34° ± 6° in one year, respectively. For both solutions, P Q , will approach its maximum 
extent, 3600, in less than 2 years. Pulsar A is expected to disappear in about 14 years 
according to solution 1. For solution 2, the pulsar will disappear in approximately 4.5 
years and then reappear in 10 years as a single-peaked pulsar (temporarily) where it 
will remain "on" for 6-7 years. For both solutions, we estimate that the pulsar came 
into view approximately 4-5 years in the past. This can explain why the Parkes 70-cm 
survey for pulsars, completed in 1997, did not detect this system in spite of sufficient 
sensitivity and appropriate sky coverage 10 . 

Recently, Demorest et al. 11 used the polarization properties of A together with the 
standard rotating vector model of Radhakrishnan & Cooke 12 in an attempt to measure 
its emission geometry. Due to limitations in both the model and the data, they were 
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only able to measure the angle a which they found to be 4° ± 3°. This result is 
consistent with our solution 1. Since they only measured one of the five angles needed 
to completely specify the geometry, they were unable to predict the future evolution of 
the emission properties of this system. 

We note that the current form of the model does not consider the precession of B. 
However, it is likely that wind-torques from A, which dominates the system energeti- 
cally, have caused the spin axis of B to align with the direction of the orbital angular 
momentum 13,14 . In this case, geodetic precession will have no effect on the emission 
from B. It is also possible that the direction of B's emission beam may lie in the or- 
bital plane as a result of the stimulated emission process regardless of the magnetic 
field alignment. In this scenario, geodetic precession would also have little effect on the 
direction of B's emission. 

Given that the model presented here accurately describes the current data, it be- 
comes important to understand the physics behind the stimulation process. Future 
work will explore various possible mechanisms. One idea involves "jump-starting" the 
pulsar emission processes in B by initiating electron-positron pair cascades in its mag- 
netosphere that emit coherent radio emission. The initiating particles could be the 
positrons and electrons emitted in the wind of A, or, more likely, gamma rays that are 
expected to be travelling in nearly the same direction as the radio photons. Alterna- 
tively, pressure from a conal A-wind could distort B's magnetosphere 14 and push its 
beam more directly into our line of sight. 
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Table 1: Measured and Inferred Properties of J0737-3039 



Measured Parameters 



Orbital Inclination 



Profile outer width 



Profile inner width 



Orbital outer width 



Orbital inner width 



Orbital Midpoint 





r , 



88.5° ± 1.5° or 91.5° ± 1.5° 

251° ± 10° 

109° ± 10° 

110° ± 10° 

22° ± 10° 

246° ± 10° 



Inferred Parameters for i = 88.5 ± 1.5 



Spin Axis 

Magnetic field angle 
Emission ring opening angle 
Emission ring width 
Total visible time 
Time remaining 

Predicted changes in 1 year 

Change in P Q 
Change in Pj 
Change in 
Change in Oj 
Change in O m 



A 

<t> 
a 

P 
25 



Solution 1 
167° ± 10° 
246° ± 5° 
1.6° ± 1.3° 
78° ± 8° 
1.9° ± 1.4° 
19 ± 2 yrs 
14 ± 2 yrs 
Solution 1 
42° ± 16° 
31° ± 8° 
0° 
0° 

4.77° ± .01° 



Solution 2 
90° ± 10° 
239° ± 2° 
14° ± 2° 
42° ± 4° 
15° ± 2° 
19 ± 2 yrs 
4.5 ± 0.1 yrs 
Solution 2 
96° ± 13° 
34° ± 6° 
0° 
0° 

4.77° ± .01° 
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Note - The measured and inferred geometric parameters for PSR J0737-3039. Timing 2 
and scintillation 4 observations have determined i to be either 88.5° ±1.5° or 91.5° ±1.5°. Two 
values are allowed since current data are unable to distinguish between i and 180° — i, although 
future timing observations should resolve this degeneracy. If a solution to the constraint 
equations exists for a given inclination, i, then a solution will also exist for an inclination of 
180° — i. The only difference between the two solutions will be that A will become 180° — A. 
In the text, the results are presented only for the case i < 90° with the understanding that 
corresponding solutions also exist if i > 90°. Note that the evolution of a given solution does 
not change when A goes to 180° — A. We estimated the errors in the inferred parameters by 
allowing the measured parameters (i, P Q , Pi, Q , Oi, O m ) to vary by the given uncertainties 
and then calculating the variations in the resulting solutions. The solutions were found by 
solving the five non-linear constraint equations using Broyden's method 15 . An initial guess 
for each of the five angles was chosen at random and then refined using this algorithm. This 
was performed 10000 times in order to determine all possible solutions. 
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Figure 1 - The hollow cone model of radio pulsar emission. 17 is the spin axis of the pulsar, 
which is separated by the magnetic dipole axis \i by the angle a. The conal emission centered 
on p has an opening half-angle p and is of angular thickness 25. The cut of an observer's 
line-of-sight through the cone as it rotates can produce either a single or double-peaked pulse 
profile. This emission geometry together with the stimulated emission hypothesis yield the 
following five constraint equations: 



sin(A) cos 



Oi 
2 



cos(p + 5 — a) ;a>p + 5 

sin(A) ;p-8<a<p + S (!) 

cos(/9 — 5 — a) ; q < p — 5 



sin(A) cos (^Y^j = cos (p + 5 + a) (2) 

(p = O m (3) 

cos(p-S) = A cos (a) + Vl- A 2 sin(a) cos (.Pi/2) (4) 

cos(/9 + (5) = A cos(a) + \/I - A 2 sin(a) cos(P D /2) (5) 

Here, A = cos(A)cos(i) — sin(z) sin(A) sin(</>). Pi, P Q , Oi, O a , O m , and i are measured 
parameters defined in the text. A and <fi are defined in figure 2. 

Figure 2 - The geometry of key orbital- and spin-related angles in the J0737-3039 system. 
Orthogonal axes representing the three spatial dimensions are denoted by x (the direction 
of the Line of Nodes in the orbital plane), y (the projected direction of the Line of Sight of 
an Earth-bound observer onto the orbital plane), and z (the direction anti-aligned with the 
angular momentum vector of the orbit, J). The pulsars orbit in the x-y plane in a clockwise 
direction. The orientation of the spin axis of J0373-3039A is 0, which is separated by the 
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orbital angular momentum, J, by the angle A, and the magnetic dipole axis, fx, by the angle 
a. The projection of f2 onto the orbital plane defines the angle 4>. The line-of-sight of the 
observer is inclined by the inclination angle i from J in the y-z plane. 

Figure 3 - A greyscale plot of the A pulsar pulse profile evolution as a function of time for 
both sets of solutions. The profile is given by a constant time slice. T=0 refers to January 1, 
2004 (MJD= 53005.0). Since this model only determines the locations of the emission regions 
as opposed to the exact shape of the profile, the simulated pulse profiles are represented as 
square waves. We predict that the 'A' pulsar will disappear from view in about 14 years for 
solution 1 or about 4.5 years for solution 2. The lack of pulsed emission from either solution 
can explain the lack of detection of the pulsars during the Parkes 70-cm pulsar survey in 
the mid-1990s 6 . Note that we have only considered one emission cone emanating from one 
side of the A pulsar. Given the symmetries of a dipole magnetic field, it is possible that 
another emission cone, the "auxiliary cone", is emanating from the other magnetic pole. If 
it exists, this auxiliary cone may eventually appear as a new component in the pulsar profile 
of A due to geodetic precession. Its appearance, though, should not affect our predictions 
for the positioning of the components of the currently observed cone. In the framework of 
the stimulated emission hypothesis, the auxiliary cone will also induce emission from B. This 
"secondary" emission should appear at orbital phases approximately 180° away from those 
at which enhanced emission is currently seen. Intriguingly, secondary emission from B has 
been reported from these orbital phases 11 , but the enhancement factor is much smaller than 
that in the two bright regions. Absorption or scattering of the radio emission at the shock 
generated near B by the interaction between the plasma winds from stars could explain this 
effect. A similar shock is believed to cause the eclipses seen in A 13 . 




Figure 1 
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